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Abstract. An irreducible algebraic decomposition Lif =0 Xi = U^_ (U -L^Xy) 
of an affine algebraic variety X can be represented as an union of finite disjoint 
sets uf_ Wi = Uj_g(U^l_ 1 Wi i j) called numerical irreducible decomposition (cf. 
[Ill^nil-ini^IIHl^IISl-Ell-ISl,!!!). W» corresponds to a pure i-dimensional X* , 
and Wij presents an i-dimensional irreducible component Xy . Modifying this 
concepts by using partially Grobner bases, local dimension, and the "Zero 
Sum Relation" we present in this paper an implementation in SINGULAR to 
compute the numerical irreducible decomposition. We will give some examples 
and timings, which show that the modified algorithms are more efficient if 
the number of variables is not too large. For a large number of variables 
Bertini is more efficient. Note that each step of the numerical decomposition 
is parallelizable. For our comparisons we did not use the parallel version of 
Bertini . 

keyword: Witness point set, Homotopy function, Grobner basis, Local dimen- 
sion, Monodromy action, Zero Sum Relation. 

1. Introduction 

Given a system of n polynomials in C^, 

f(x 1 ,...,x N ) := 

Let X=V(f) be the algebraic variety defined by the system above. X has an 
unique algebraic decomposition into d pure i-dimensional components Xi, X = 
uf =0 AV Where X t = U^A%- is the union of di i-dimensional irreducible compo- 
nents, do, di, dd positive integers. 

The numerical irreducible decomposition (cf. [15| , [17] , [18] , [19] , [22] ) is given as the 
union W = Uf =0 Wi = uf =0 (\jf =1 W^). The W t are called i- Witness point sets and 
are given as an intersection of the pure i-dimensional component Xi of X with a 
generic linear space L in C N of dimension N-i, the Wij are called the irreducible 
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witness point sets presenting the irreducible components Xij of dimension i with 
the following properties: 

• Wij consists of a finite number of points contained in . 

• tt(Wij) =deg(X ij ) fori^O. 

• W i:j n W a = for j =/= I. 

Computing of the numerical irreducible decomposition uses numerical polynomial 
homotopy continuation methods(cf. 24 , 25 ). This requires that the number n 
of a given polynomial system has to be equal to the number N of the variables. 
Therefore we reduce the polynomial system which defines X to a square system of N 
polynomials in N variables (cf. [15] , [17] , [22] ) . Numerical irreducible decomposition 
(cf. [T5],[T7],[I5],[T5],[22]) is proceeded in three steps: 

• I s * step reduces the polynomial system to a system of N polynomials in N 
variables and computes a finite set Wi called witness point super set for a 
pure i-dimensional component Xi for i=0,...,N-l. Wi consists of points on 
Xi and Ji a set of points on components of larger dimension the so-called 
Junk point set (cf. [15 ] . [2"2 ] ) . 

• 2 nd step removes the points of Ji from Wi to obtain a subset Wi of the pure 
i-dimensional component Xi (cf. [2"2"]V 

• 3 rd step breakup Wi into irreducible witness point sets representing the 
i-dimensional irreducible components of X using two algorithms. The first 
algorithm finds points on the same irreducible component in the Witness 
point set connected by path tracking techniques applying the idea of mon- 
odromy. The second algorithm computes a linear trace for each component, 
which certifies the decomposition (cf. |14 | .|21 | ). 

In the second section we present a modified algorithm to compute Wi using the 
cascade algorithm (cf. [15] , [17] , [22] ) and Grobner bases in the zero-dimensional 
case. Then we use the homotopy function (cf. [H]), local dimension and Grobner 
bases in the zero-dimensional case to remove Junk points from Wi to obtain the 
i- Witness point set Wi. This is explained in third section. In the fourth section 
we explain how to use the "Zero Sum Relation" and monodromy action on the 
algebraic variety to breakup Wi into irreducible witness point sets. In the fifth sec- 
tion some examples are tested and timings are given on the basis of our Singular 
implementation. 

2. Witness Point Super Set 

Definition 2.1. Let Z=V(f) be an affine algebraic variety in C N of dimension d, 
and X be a pure i-dimensional component of Z. Let Li be a generic linear space in 
C N of dimension N-i. A set Wi C is called i- Witness Point Super Set for X if 
it has the following properties: 

• Wi is a finite set of points. 

• x nLi cf, c z n Li . 

The union W of all i- Witness Point Super Sets is called a Witness Point Super Set 
for Z. 

In (cf. [TS] , [T7] , [32] ) the cascade algorithm is used to compute Wi. It starts 
with i=N-l to compute the Witness point Super Sets Wi- It needs to define a start 
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system G(x)=0 for the homotopy continuation method (cf. [H],[1S]) and needs to 
know its solutions. We use a Grobner basis to compute the dimension d of Z, then 
use the cascade algorithm (cf. [TS] , [T7] , [32] ) which starts with i=d-l. We show 
below that we do not need to define a start system. 

Algorithm 1 WitnessPointSuperSet 
Input: F x , ...,F n <E C[xi, ...,x N ]. 

Output: .., f N }, {W r , .., W d },L. {f lt .., f N } square system, Wj Witness point 
super sets corresponding to a pure i-dimensional component of V(/i, ...,Jn), L 
a set of linear polynomials defining a linear space of dimension N-d. 

/ = {/i,-..,/jv} reduction of F = {Fi,...,F n } to a square system 

(<*mmmy, 

d=dim(y(/i, /jv)), using Grobner basis (cf . [8] . [TT] ) ; 

r=N-rank(f jj, rank(f) the rank of the Jacobian matrix of the system f at a generic 
point; 

L = {li, ...,ld} a set of d generic linear polynomials; 
if d = r then 

compute Td = V{f\, /jy, h, ■■■> Id), using triangular sets, (cf.[5],[TT]); 
set W d = {{xi,...,x N ) \Jxi,...,x N ) £ T d , (xi,...,x N ) € V(F)}; 
return {f x , f N }, {W d }, L ; 
else 

for i — r to d do 
if i = then 

else 



Cli(f)(x,Zi,...,Zi) =: 



fN{x) +J2] = l X N J Z j 

h + z\ 



V h + Zi I 

Xkj 6 C generic, k = 1, N, j = 1, i; 
for i — d to r do 
if i = d then 

compute Tj = V(fij(/)(a;, zi, 2^)), using triangular sets, (cf. [8].[TT]); 
set Wi = {(xi,.. .,Xff) | (xi,...,x N ,0,...,0) eTi, (xi,...,x N ) eV(F)}; 
S i = T i \{(xi,...,x N ,zi,...,z i )eT i | z\ = .... = Zi = 0}; 
else 

compute Tj = V(Qi(/)(x, zi, z^)), using homotopy function with 
Qi+i(f)(x, z%, Zi) as start system and Si+i as start solution set (cf. 
BSI.IIZIM; _ _ 
return {fi,...J N }, {W r , ...,W d }, L ; 



V(fi, /jv) has no components of dimension smaller then N — rank(f) (cf. [22]). 
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3. Computation Witness Point Set 

The witness point super set Wi is an union of an i- Witness point set Wi and a 
junk point set J* (cf. [T5] , [T7] . [22] ) . 

Wi = W. l U Ji, where Wi C X t and J l C Uj >t Xj for i = 0, 1, d. 

We use Grobner bases for the 0-dimcnsional ideal, local dimension and homotopy 
continuation method to remove the points of Ji from W, as follows: 



Algorithm 2 WitnessPointSet 

Input: {/i,-.,/iv} C C[xi, .., xn], {W r , .., Wd} a list of witness point superset sets, 
L = .., Id} a set of generic linear polynomials (Output of Algorithm 1). 

Output: {/i, /jv}, {W r , ■-, Wd}, L = {li, Id}- Wi a Witness point set corre- 
sponding to a pure i-dimensional component of V(fi, /at). 

W d = W d , s d = Wd\ 
for i — d — 1 to r do 

W t = W; 

for each point w £ Wi do 

compute^ t = dim w Z for Z = V(f\ — ,fi(w), fx — fjsr(w)), using Grobner 

basis (cf. [H],[n]); 

if t > i then 
W t = W % \ {w}; 
for each point w £ Wi do 

if i = then 

choose A C C dxN a generic matrix and a generic e G C , ||e|| small; 
compute S = V({fi, Jn, A(x - w)}), T = V({fi, Jn, A(x -w - e)}), 
using triangular sets (cf. [8].[TT]): 
if $S = jJT then 

Wi = Wi \ W; 

else 

for j = i + 1 to d do 

choose A C xN a generic matrix; 
if j = d then 

compute S = V({f±, /at, A(x — w)}), using triangular sets, (cf. 

0.E1); 

if jjS 1 = Sd then 
= Wi \ W; 

else 

compute S — V({fi, Jn, A(x — w)}), using homotopy function with 
start system {fx, /jy, h, lj} and start solution Wj (cf. [23]); 
if w £ S then 
W l = Wi \ {w}; 
return {f x , f N }, {W r , W d }, L ; 



2 t = dim w Z is less or equal to the local dimension of a point v g V(/i, ...,/jv) (cf. [9]) with 
\\w — v\\ small. 



4. Partition Witness Point Sets 



In this section we show that the monodromy action on an algebraic variety Z 
and the Zero Sum Relation are sufficient to find the breakup of the k- Witness point 
set Wk into irreducible k- Witness point sets. We present here a modified version of 
the algorithms described in (cf. [14].|21]). 

Let Z be a pure k-dimensional algebraic variety in C N , and Z — \J r i=1 Zi be the 
irreducible decomposition of Z. Let n : C N — > C k be a generic projection and let 
/ C C fc be a general line. 
Set 

• W/ := 7T — 1 (Z) n Z a set of r different curves in C^. 

• U is a non-empty open subset of I consisting of all points x £ I with t:^ 1 (x) 
transverse to Z. 

• W := tt^ 1 (x) H Z for a generic element x £ U, and V x := V a subset of W. 

• Wi :— tt^ 1 (x) n Zi for an irreducible k-dimensional component Zi of Z. 

• A : C N — > C a linear function one-to-one on W . 
Define the function s : U — > C by 

s(y) = A 0)- 

ZGVy 

Where V y is a subset of n~ 1 (y) n Z defined by. 

V y := {z | z on a curve through a point of V}. 




x y 



Theorem 4.1. Let I, U, W, Wi for i=l,...,r, and the functions X, s be as above. 
If the function s is continuous and VnWi ^ for some i £ {1, r} 7 then Wi C V. 

Example 4.1. Before proving the theorem we illustrate it by an example. 

• Let Z be the algebraic variety of dimension one in C 2 defined by the poly- 
nomial f(x,y) = (x 2 + y 2 — 5)(x — 2y — 3). Let L\ be the linear space of 
dimension one in C 2 defined by the polynomial l\ — x + y — 3. 

• Define a homotopy function : 

a : [0,1] — >p~ l {L) given by 

a(t) = (1 - t)l +th=x + y-2t-l. 

Let Lq be the 1-dimensional linear space defined by the polynomial Iq = 
x + y — 1, Lq n Z 7^ 0. Let G(N — k, N) be the Grassmanian and R := 
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{(L N - k ,x) G G(N - k,N) x ~Z | x G Ljv-fc H Z} fee t/ie family of the 
intersections L_/v_fc PI Z , Ljy-k C P w k- dimensional linear spaces and Z C 
ifte closure of Z. Letp : R — >■ G(N — k, N) be the canonical projection. 

Then with conditions above a(t) maps a point in L±n Z to a point in Lq l~l Z as t 
goes from 1 to 0. 

Proof, (of theorem 4.1) Assume that Wi V . Since Wi DV ^ 0, then there are 
a,b 6 such that a is not in V and b G V. Let ai,...,a r denote the points of 
the set V \ {&}. By (Corollary 3.5 in [14]) there is a loop a in the fundamental 
group tvi(U, tt~ 1 {x)) with a(0) = a(l) which takes dj to <Zj for all j=l,...,r, and 
interchanges a and b. 

Since a is a continuous loop and s : U — 
s o a : [0,1] — > C is continuous and 



is continuous, the composition 



a(a(l)) = s(a(0)) 



A(a) 



A(6) 



as t goes from 1 to 0. This implies that A (a) 
that A is one-to-one on W. Thus Wi C V. 



A(6). But this contradicts the fact 

□ 



Example 4.2. Let Z be a pure 1-dimensional component in C 2 defined be the poly- 
nomial f(x, y) = (y — x)(y — 2x)(y — 3x), and Z = Z\ U Z2 U Zj, be the irreducible 
decomposition. Let it : C 2 — 5- C be the projection given by n(x, y) = x, and 
A:C 2 — >C, X(x,y) =y. 

Note that the restriction ofir to Z, itz is proper and generically three-to-one with de- 
gree 3 equal to the degree of Z. A is one-to-one on the fiber it" 1 {y) — {(x,x), (x,2x), 
(x, 3x)}. Let L be the linear space in the Grassmannian G(l,2) defined by the 
linear polynomial l(x,y) — x + y — 2. L intersects Z in the finite set W := 
{(1,1), (f, |), (|, I)}- 

LetV := {(1,1),(|,|)} C W. The function J2 ve v X ( v ) 

given by X(x, x)+X(x, 2x) = 

x + 2x = 3x is continuous. Then by the theorem above for an irreducible 1-Witness 
point set W\ CV it contains {(1, 1)}. 

Preparation of the algorithm(IrrWitnessPointSet): 

Let Zu = U r i=1 Zki be the union of the irreducible k-dimensional components of the 
algebraic variety Z and Lk be the generic linear space in defined by k linear 
equations 

lj — Cj + Cj\X\ + ... + CjNXN- 

for j=l,..,k and i=0,l,...,N , dj E C. 

We use the generic linear space Lk to define the generic projection 



Z k+1 , w(xx,...,x N ) := (z 1 ,...,z k ,z k+1 ) as follows: 



\ x N j 



( 



I— > 



\ 



V z k +i J 



( Cn C12 

C21 C22 



Cfcl C k 2 

V Pi P2 
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ClN \ ( Xl \ 

ClN 



CkN 
PN ) 



\ X N J 



Pi, ...,Pn £ C randomly chosen. 

Set X(xi, ...,xn) := Zk+i and I :— z k C C as in the theorem above. Let y := Cko 
vary and fix the other cio, Cfk-i)o- 

• Let Lfc.j, be the linear spaces defined by the linear equations l\,...,lk-\ 
above and Zfc i3/ := y + c k \X\ + ... + c kN x N . 

• For the subset V y of W y := L^ y n Zfc the function s is given as a function 
s y : Z k C\ Lk, y — > C, 

s y (xi, ...,£Civ) := A(xi, ...,xat). 

(ssi,...,a:jr)eV5/ 

It is convenient to test the linearity of s y . 

• Since y = —{ck\X\ + ... + c^xm), then we can define a function s 

s:C — >• C, s(y) := s y (x\,...,x N ) = ^ A(:Cl, ...,a:jv)■ 
(:^:l,...,:^;iv)6^ / H 

• s y is a linear function in x = (x\, 2 at) if and only if s is linear in y. 

• To test the linearity of s, we take three values of y in C, say a, b, c. 
If there exist A, B £ C such that: 

(s(a) = Aa + B, s(b) = Ab + B) => s(c) = Ac + B. (4.1) 

Then s is linear in y. 

Here s(a), s(b) and s(c) correspond to the subsets V a C W a = Z k n Lk, a , 
Vb C W& = Zfc n Lfc^ and V c C W c = Zk n Lfe iC respectively with ftV a = 
ttH = Wc = m. 

So far this is the approach which can be found in [14J . Now we give some modifi- 
cations. 

• The condition (4.1) of the linearity above is equivalent to the following 
equation 

s(a)(b -c) + s(b)(c -a) + s(c)(o - b) = 0. (4.2) 

• From Theorem 4. 1 we obtain: If Wkj fl V a ^ for some j £ {1, r} and 
the condition (4.1) of the linearity above is true, then Wkj Q V a - 

• Let Z(y) := {z = Y.t=iPt v t I v = (vi, ...,v N ) £ V y ,p = (p x , ...,p N ) £ C N }. 
Then 

N 

s(y)= = £ = £ *. 

• To compute the sets VJ, and 14 we use the homotopy function as t goes from 
1 to using V a as start set. In particular 

V b := ((1 - t)L Kb + tL k , a ) n Z, V c := ((1 - t)L k , c + tL k , a ) n as t goes 1 -> 0. 

Continuation of the homotopy function implies that the i-th points in the 
sets V a , Vb and V c are on the same irreducible component. 

• Let V a := {vi,...,v m }, V b := and V c := {vi,...,v m } be the sets 
computed by using the homotopy function above . Let Z(a) := {a%, a m }, 
Zip) := {&i,..., 6 m } and Z(c) := {ci,...,c m } be the sets corresponding to 
the set V a , Vb and V c respectively and defined as Z(y) above. 
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From (4.2) we get a condition equivalent to the condition (4.1) of the lin- 
earity 

tci rn m 

(b-c)J2 <H + (c-a)J2 b i + ( a - b )J2 °i = °- ( 4 - 3 ) 

i—1 i—1 i—1 

The condition (4.3) is called Zero Sum Relation (cf. [7]) of a given subset 
V a QW denoted by ZSR(V a ). 
• The sets V a , V& and V c have distinct points and the same cardinality m, 
then obviously 

ZSR(V a ) = ]T ZSR({ ai }). (4.4) 

ffl,£V 

where ZSR({a,i}) = (b — c)aj + (c — a)bi + (a — b)ci is defined as a Zero 
Sum Relation of a given point in V a . 

Algorithm 3 IrrWitnessPointSet 

Input: {/i, /jv} C C[xi, ...,xjv], {W r , Wd}, a list of witness point sets, L — 
{h,.-.,ld} a set of generic linear polynomials (Output of Algorithm 2). Where 
Wk = {wi, ...,w mk } are witness point sets for a pure /c-dimensional component 
Z k ofZ = V(f 1 ,...,f N ),k = r,...,d. 

Output: {{W r i,...,W rtr }, {Wai,...,Wdt d }}, W krk irreducible Witness point 
sets corresponding to a fc-dimensional irreducible component Z krk of Z^ . 

for k = r to d do 

a:=CfeOj define Lk a to be the linear space defined by the subset {l±, Ik} C L. 

choose b, c G C generic, define L^, Lfe c as above; 
I^ a = ^ fe , W^ fc = 0, W c = 0, R = $; 
choose pi, ...,pn 6 C; 
for i = 1 to m/j do 

compute {vi} C Z n Lj^j, and {di} C2fl Lfc ;C using the homotopy function 
with {/i, /jy, ?Xj •••) ife— i) ^fe,a} as start system and {wi} as start solution; 
compute the Zero Sum Relation of {w;}: 

N N N 

n=(a- b)(Yp 3 v l3 ) + (b- cj^pfWij) + (c - a)(^p t ^); 
i=i j=i i=i 

fl= J RU{r i ]0; 
int 4^ = 0; 
while R ^ do 

if X)teT ^ — an d T is a smallest subselQof i? then 

tfc=ik + l; 

M4t fc C W consists of the points corresponding of the points of T; 
R = R\T; 

return {{W rl , W rtr }, {W dl , W di<i }} ; 



the i — th point in R corresponds to the i — th point in W a ; 
4 smallest subset with respect to the cardinality. 
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We give an example of a pure 2-dimensional variety Z which is an union of two 
2-dimensional irreducible components Z\ and Z2- Z\ is of degree three and Z<i 
is of degree two. The 2- Witness point set W for Z is given as a finite subset of 
Z consisting of five points {u>i, W2, W3, 1114, w§}. Z\ should contain three points 
W\ := {wi, U)2, W3} and the remaining points W2 '■= {w^Wq} are on Z2. The 
algorithms (cf. [14],[2T]) use the homotopy function at least nine times to breakup 
W into W\ and W2 ■ We will show below that we do not need more than five times 
to use the homotopy function to breakup W into W\ and W2 ■ 

Example 4.3. . 

Let Z be the algebraic variety of dimension two in C 3 defined by the polynomial 
f(x,y,z) — (x 3 + z)(x 2 — y). Let L be the linear space of dimension one in C 3 
defined by the linear equations l\ = Ax + 7y + 2z + 6, I2 — 5x + 7y + 3z + 6. Then 
W := LDZ = {wi,w 2 ,W3,Wi,w 5 }, wher$ Wl = (1,-1.1428571429,-1) 
(0,-0.8571428571,0), 

w 3 = (-0.1428571429 + i * 0.9147320339, -0.8163265306 - % * 0.2613520097, 
0.1428571429 - i * 0.9147320339), 
w 4 = (-1, -0.5714285714, 1), 

w 5 = (-0.1428571429 - i * 0.9147320339, -0.8163265306 + i * 0.2613520097, 
0.1428571429 + i * 0.9147320339) . 

We now illustrate Algorithm3 (IrrWitnessPointSet) : 

• Use the linear space L\ to define the linear projection ir : C 3 — > 
follows 



,w 2 




w(x,y,z):=\ 5 7 3 \ \ y \= (4a; + 7y + 2z, 5x + 7y + 3z,x + 2y + 3z). 



Define the linear space Li iC of dimension one in C 3 by the linear equations 
l\ — 4x + 7y + 2z + 6, l c — 5x + 7y + 3z + c, where c is generically chosen 
in C. Then 

nznL 1:C ( x i V) = (—6, ~c, x + 2y + 3z). 
Define the linear function A : C 2 — > C by X(x, y, z) :— x + 2y + 'iz. 
For a=6, let Vi = V a := {w n = (1,-1.1428571429,-1)} C W, L ha := 
L the linear space defined by li = 4x + 7y + 2z + 6, l a = hx + 7y + 
3z + 6. Then5Z(a) = {J2 ve v a x (v) = ^nW + 2(t»n[2]) + 3(tfl U [3])} = 
{-4.2857142858}. 

Let b—9, L\j> the linear space defined by l\ = Ax + 7y + 2z + 6, If, — 
5x + 7y + 3z + 9. Compute V b := {tL lM + (1 - t)L ltb ) n Z = {w 12 = 
(1.671699881657157, -0.4776285376163331, -4.671699881657164)} as t goes 
from 1 to 0, using V a as a start solution. Z(b) — {1012 [1] + 2(u>i2[2]) + 
3(wi2 [3])} = {-13.2986568385470012}. 

Let c=63, L\ tC the linear space defined by l\ = Ax + 7y + 2z + 6, l c = 
5x + 7y + 3z + 63. Compute V c :— (tL\_ a + (1 — t)L\_ c ) D Z = {wis — 
(3.935100643260828, 14.30425695906836, -60.93510064326094)} as t goes 



^Notc that the values of Wi are approximate values. The following equalities are therefore to 
interpret as approximations of the points Wi. 

"we use the notation Wij = (viij -[1], Wij [2], Wij [3]) for i=l,..,5, j=l,2,3. 
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from 1 to 0, using V a as a start solution. Z(c) = {wi3[l] + 2(wi 3 [2]) + 
3(wi 3 [3])} = {-150.261687368385272}. 

Zero Sum Relation of V x = {(1, -1.1428571429, -1)}; 

n : = e ( b c ) + E ( c - «) + E ( a - 6 ) = 

aGZ(a) beZ(b) ceZ(c) 

= -75.8098062588232524. 

The zero sum relation set of Vi = {(1,-1.1428571429,-1)} is R x := {n = 
-75.8098062588232524}. 

• Let a=6, V a := {w n = (0,-0.8571428571,0)} C W, L 1>a := L the linear 
space defined by l\ = 4x + 7y + 2z + 6, l a = 5x + 7y + 3z + 6. Then 
Z{a) = {E ve v a A («) - wii[1]+2(«;ii[2]) + 3(i«ii[3])} ='{-1-7142857142}. 

• Let b=9, L\fi the linear space defined by li = Ax + 7y + 2z + 6, h = 
5x + 7y + 3z + 9. Compute V b := (tL lM + (1 - t)L hb ) n Z = {w 12 = 
(-0.8358499408285809 + i * 1.046869318849985, 0.2388142688081706 - i * 
0.2991055196714253,-2.164150059171436 - i * 1.046869318849981)} as t 
goes from 1 to 0, using V a as a start solution. Z(b) — {wi 2 [l] + 2(^12 [2]) + 
3(wi2 [3])} = {-6.8506715807265477- i * 2.6919496770428086}. 

• Let c=6S, L\ tC the linear space defined by l\ = Ax + 7y + 2z + 6, l c = 
hx + 7y + 3z + 63. Compute V c := (tL lM + (1 - t)L l>c ) n Z = 

{w 13 = (-1.967550321630417+ i * 3.257877039491183, 
15.99072866332302 - i * 0.9308220112831772, 

-55.03244967836969- i * 3.257877039491242); } as t goes from 1 to 0, us- 
ing V a as a start solution. Z(c) = {w\ 3 [l} + 2(wi 3 [2]) + 3(wi 3 [3])} = 
{-135.083442030093447 - i * 8.3773981015488974}. 

Zero Sum Relation of V 2 = {(0, -0.8571428571,0)}: 

r 2 := E ( & - c )+ E ( c -«)+ E ( fl - & ) = 

aeZ(a) beZ(b) c£Z(c) 

= 107.3334745556671221 - i * 128.308937286793398. 

The zero sum relation set of V 2 = {(0, -0.8571428571,0)} is 
R 2 := {r 2 = 107.3334745556671221 - i * 128.308937286793398}. 

• For the other points V 3 = {u> 3 }, V4 = {W4} and V5 = {W5}, we found the 
zero sum relations R 3 := {r 3 = -9.38237104997583366+i*127.0170767088} ; 
R 4 := {r 4 = -31.5236682999307779+ i * 128.3089372867945956} and 

R 5 := {r 5 = 9.382371038077068 - % * 127.0170767088}. 

• The set of Zero Sum Relation for all points of W is R = Li j =1 Rj = 
{ri, r 2 , r 3 , r§}, where i-th point in W corresponds i-th point in R. 

• Find the smallest subset T of R with J2teT t = 0, which corresponds an irre- 
ducible Witness point set of W. Then we get T\ — {r 3 , r 5 }, T 2 = {n, r 2 , r 4 } 
corresponding to the irreducible Witness point sets W\ — {w 3l w§}, W 2 = 
{wi,w 2 ,u>4} respectively. 

Remark 4.1. The points of a Witness point set are computed approximately by 
using the homotopy continuation method. Therefore the result of the Zero Sum 
Relation is only almost zero. 
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5. Examples and timings with Singular and Bertini 



In this section we provide examples with timings of the algorithms 
WitnessPointSuperSet, WitnessPointSet, and IrrWitnessPointSet implemented 
in Singular to compute the numerical decomposition of a given algebraic variety 
defined by a polynomial system and compare them with the results of Bertini . 
We did not use the parallel features of Bertini . 

We tested to versions of the implementations in Bertini using the cascade al- 
gorithm and using the regenerative cascade algorithm. Timings are conducted by 
using the 32-bit version of Singular 3-1-1 (cf. |S]) and Bertini 1.2 (cf. 3 ) on an 
Intel® Core(TM)2 Duo CPU P8400 @ 2.26 GHz 2.27 GHz, 4 GB RAM under 
the Kubuntu Linux operating system. 



Let Z be the algebraic variety defined by the following polynomial system: 
Example 5.1. (cf. [17] ). 

(y - x 2 )(x 2 + y 2 + z 2 - l)(x - |) 
(z - x 3 )\x 2 +y 2 + z 2 - l)(y - |) 
(y - x 2 )(z - x 3 )(x 2 + y 2 + z 2 - l)(z - §) 

Example 5.2. (cf. [22], Example 13.6.4). 



f(rvz)-l x(y 2 - x 3 )(x - 1) 
T[x,y,z)-\ x ( y 2 ^ x 3^ y ^ 2 ){3x + y) 

(x 3 + z)(x 2 — y) 
f(x,y,z) = \ (x 3 +y){x 2 - z) 

{x 3 + z)(x 3 + y)(z 2 — y) 

x(y 2 — x 3 )(x — 1) 
f(x, y,z)= \ x(3x + y)(y 2 - x 3 )(y - 2) 
x(y 2 — x 3 )(x 2 — y) 

(x - l)((x 3 + z) + (x 2 - y)) 
f(x,y,z)= \ (x 3 + z)(x 2 ~y) 

(x 3 + z){x 2 - 1) 

Example 5.6. 

/ (y - x 2 )(x 2 + y 2 + z 2 - l)(x - |) + x 5 
f(x, y,z)=[ {z- x 3 )(x 2 +y 2 + z 2 - \){y - |) + y A 

\{y- x 2 )(z - x 3 )(x 2 +y 2 + z 2 - l)(z - |) + z e 

Example 5.7. 

x(y 2 - x 3 )(x - 1) + y 2 
x(y 2 - x 3 )(y - 2)(3x + y) + x 



Example 5.3. 



Example 5.4. 



Example 5.5. 



f(x 1 y,z)— ( , 2 ; . .),,•.,. ... 



li 



Example 5.8. 



Example 5.9. 

( 



f{x,y,z) = 



(x 3 + z){x 2 -y) + x 4 
(x 3 + y)(x 2 - z) + y 3 
{x 3 + z){x 3 + y)(z 2 - y) + 



f{x,y) = 



/i = -3568891411860300072a; 5 + 1948764938x 4 + \ 
3568891411860300072xV - 1948764938a;?/ 2 



f(x,y,z) 



f(x,y,z) 



+ 2x 3 y 2 z 2 + 2xy 4 z 2 + xy 2 z 4 



-xy' 

: - 2x 3 y 2 z - xy 4 z - 2x 2 y 2 z 2 - 2y 4 z 2 



A 2 

- x y 



z 5 - 3x 3 y 2 - 3xy 4 + x 4 z+ 



Example 5.11. 

/ h = x 5 y 2 + 2x 3 y 4 - 
-2x 2 y 4 -y 6 - x 5 z 
2x 3 z 3 - 2xy 2 z 3 - y 2 z 4 - xz 
2x 2 y 2 z + y 4 z - 3xy 2 z 2 + 2x 2 z 3 + 2y 2 z 3 + z 5 + 3x 2 y 2 -r 
3?y 4 + 2>x 3 z + 2>xy 2 z + 2>y 2 z 2 + 3xz 3 + 2xy 2 - 3x 2 z - 3y 2 z- 
3z 3 - 2y 2 - 2xz + 2z 

h = x e y + 2x 4 y 3 + x 2 y 5 + 2x 4 yz 2 + 2x 2 y 3 z 2 + x 2 yz 4 - 



5x b - 10xV 



I0x 2 y 2 z 2 



— 5x 2 y 4 — x 4 yz - 
2x 2 yz 3 - 2y 3 z 3 



2x z y 6 z 
h 4 



5x 2 ~ 4 



y 5 z 



Wx 4 



yz b — 3a; y- 



3x 2 y 3 + 5x 4 z + 10x 2 y 2 z + 5y 4 z - 3x 2 yz 2 + 10x 2 z 3 + 
10y 2 z 3 + 5z 5 + 15a; 4 + 15x 2 y 2 + 3x 2 yz + 3y 3 z + \hx 2 z 2 
+3yz 3 + 2x 2 y - 15x 2 z - 15y 2 z - 15z 3 - 10a; 2 - 2yz + lOz 

h = x 6 y 2 z + 2x 4 y 4 z + x 2 y 6 z + 2x 4 y 2 z 3 + 2x 2 y 4 z 3 + 
x 2 y 2 z 5 - 7x 6 y 2 - Ux 4 y 4 - 7x 2 y e - x 6 z 2 - \7x 4 y 2 z 2 - 



17x 2 y 



A z 2 



y 6 z 2 



2x 4 z 4 - \\x 2 y 2 z 4 - 2y 4 z 4 



„2„6 



y 2 z 6 + 7x 6 z + 18x 4 y 2 z + 18x 2 y 4 z + 7y 6 z + 15xV 3 + 
27x 2 y 2 z 3 + \hy 4 z 3 + 9x 2 z 5 + 9y 2 z 5 + z 7 + 21x 4 y 2 + 
21x 2 y 4 - 4x 4 z 2 + 13x 2 y 2 z 2 - 4y 4 z 2 - llx 2 z 4 - llyV- 
7z 6 - 21x 4 z - 40x 2 y 2 z 01 - 4 ~ 0^2,3 o/i„.2„3 q,5 



21y 4 z 



24x 2 z 3 



24?y 2 z 3 - 3z 5 



14x 2 y 2 + 19x 2 z 2 + 19y 2 z 2 + 21z 4 + 14x 2 z + 14y 2 z+ 



2z 3 - 14z 2 



f 2 = -5105200242937540320x 5 ?y - 1701733414312513440xV + 
11692589628x 5 + 3897529876x 4 ?y + 5105200242937540320xV + 
\ 1701733414312513440xy 4 - 1 1692589628a; V - 3897529876xy 3 J 

Example 5.10. 

( ,h = -356737285367005125x 5 - 92300457164036000x 3 y+ \ 
1121648050080163317x 2 z + 290209720279281056y2 



h = -356737285367005125x 5 + 887060318883271500x 3 z+ 
1121648050080163317x 2 y - 2789081819567309964^ 

h = -356737285367005125x 5 z 2 + 356737285367005125x 5 y+ 
887060318883271500x 3 ^ 3 - 887060318883271500x 3 yz+ 
1121648050080163317x 2 z 3 - 1121648050080163317x 2 j/z- 
\ 2789081819567309964z 4 + 2789081819567309964jyz 2 J 



\ 
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Example 5.12. 



x\ + x\ + x 2 + Xs + x 4 - x 5 
x\ + Xi + x 2 + x 3 — X4 + x 5 

x\ + X\ + X2 — X 3 + X4 + X 5 
x\ + X\ — X2 + X 3 + Xi + X 5 
x\ — X\ + X2 + X 3 + X4 + X 5 



4 \ 



f(xi,X 2 ,X 3 ,X4,,X 5 ) 



4 
4 
4 



V 



4 / 



Example 5.13. 



/ a 2 + 2de + 2c/ + 2bg + a \ 



f(a,b,c, d,e,f,g) 



2ab + e 2 + 2df + 2cg + b 
b 2 + 2ac + 2ef + 2dg + c 
2bc + 2ad + f 2 + 2eg + d 
c 2 + 2bd + 2ae + 2fg + e 
2cd + 2be + 2a f + g 2 + f 



\ d 2 + 2ce + 2b f + 2ag + g J 
Example 5.14. . cyclic 4-roots problem. (cf\±>\,\6\). 
Example 5.15. . cyclic 5-roots problem. (cf.\5\,\§\). 
Example 5.16. . cyclic 6-roots problem. (cf.^\,^\). 
Example 5.17. . cyclic 7-roots problem. (cf.\E\,\Q\). 
Example 5.18. . cyclic 8-roots problem. (cf.\5\,\fi\). 
Example 5.19. . 

f(xu, Xl2, Xi3, Xii, Xi 5 ,X 2 l,X 2 2, X 23 ,X 2i , X 25 , ^31 , ^32 , ^33, ^34, £35) = 



TableTT] summarizes the results of the timings to compute the numerical decom- 
position.]. 



Remark 5.1. The timings show that for an increasing number of variables the 
original method of fcf. [14] . [15] . |17j . [21j . [22] ) becomes more efficient. One reason is 
that the computation of triangular sets which is used in SINGULAR for solving poly- 
nomial systems is expensive in this case. Therefore the Algorithml. Algorithm2 
become slow in this situation. This is not true for Algorithm3. 
Replacing the solving of polynomial system using triangular sets by homotopy func- 
tion methods but keeping the computation of the dimension and starting in this 
dimension is more efficient in a case of a large number of variables. 



( -X\ 2 X21 + x n x 22 \ 



-x\zx 22 + x 12 x 23 

-XUX23 + Xi 3 X 2 4 
-Xi 5 X 2 4 + X U X 25 
~X 22 X 3 i + X 21 X 32 

-x 23 x 32 + x 22 x 33 

-Z24Z33 + 3?23£34 



\ --^25^34 + £242:35 J 




(re) means using the regenerative cascade algorithm instead of the cascade algorithm 
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Example 


Bertini 


Bertini (re) 


Singular 


It i I 
10.11 


19/1 A C Q 


oys 


07 


10. Zl 


o.UoS 


9 

z.os 


1 Ada 


IQ.ol 


Imin 01 98a 

imin zi.zos 


97 /la 

z / .^s 




10.41 


lO.OOS 


9 7q 
Z. ( o 


1 77a 
I.MS 


10.01 


1 ^ QAc 
10.OOS 


o.OS 


1 9Qa 

i .zys 


10.01 


4min 13s 


15min 2s 


2min 27s 


Ic, 7I 
10. <l 


l.ooS 


1 fla 

l.DS 


u.oys 


|c 01 
IQ.ol 


3min 29s 


lOmin 43s 


i .oys 


|C Q| 

lO.al 


1 £o 

IDS 


/ S 


9a 

zs 


I0.1UI 


2niin 57s 


9Sq 

ZOS 


2min 35s 


15.111 


44min 56s 


2min 37s 


4min 3s 


15.121 


4.73s 


6s 


0.37s 


15.131 


5.84s 


8s 


Is 


15.141 


1.43s 


4.3s 


0.79s 


15.151 


3.54s 


10s 


0.57s 


15.161 


3min 23.26s 


2min 29s 


1.43s 


15.171 


2h llmin 57s 


32min 17s 


stopped after 5h 


15.181 


19h48min 17s 


6h45min2s 


stopped after 50h 


15.191 


lmin 57s 


51s 


stopped after 3h 



Table 1 . Total running times for the computing a numerical de- 
composition of the examples above 
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